********************************************************************************
*																			   *
*						Sir! I'd Rather Go to School Sir!					   *
*																			   *
*								 Mahdi Majbouri								   *
*																			   *
*									 Table 2								   *
*																			   *
*									6-6-2023								   *
*																			   *
********************************************************************************


use "census11-all.dta", clear

keep if sons==1
keep if age >= 20

keep if fage>=fage18 & fage18>=threshold - 9 & fage18<=threshold + 10


*** Regressions for sole sons
regress college_plus t fa tfa father_ed mom_ed fage mage if relate==3 & son==1 & sons == 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("SoleSons-p1") replace

regress college_plus t fa tfa fasq tfasq father_ed mom_ed fage mage if relate==3 & son==1 & sons == 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("SoleSons-p2") append

regress college_plus t fa tfa fasq tfasq facb tfacb father_ed mom_ed fage mage if relate==3 & son==1 & sons == 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("SoleSons-p3") append


*** Regressions for sole sons's sisters
regress college_plus t fa tfa father_ed mom_ed fage mage if relate==3 & son==0 & sons == 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("SoleSonSisters-p1") append

regress college_plus t fa tfa fasq tfasq father_ed mom_ed fage mage if relate==3 & son==0 & sons == 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("SoleSonSisters-p2") append

regress college_plus t fa tfa fasq tfasq facb tfacb father_ed mom_ed fage mage if relate==3 & son==0 & sons == 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("SoleSonSisters-p3") append



********************************************************************************
**  Robustness check - Regressions for boys in multiple son families
********************************************************************************

use "census11-all.dta", clear

keep if fage>=fage18 & fage18>=threshold - 9 & fage18<=threshold + 10

** Regressions for sons in multiple-son households
regress college_plus t fa tfa father_ed mom_ed fage mage if relate==3 & son==1 & sons > 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("MultipleSons-p1") append

regress college_plus t fa tfa fasq tfasq father_ed mom_ed fage mage if relate==3 & son==1 & sons > 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("MultipleSons-p2") append


regress college_plus t fa tfa fasq tfasq facb tfacb father_ed mom_ed fage mage if relate==3 & son==1 & sons > 1 & fage>=fage18 & urban == 2, vce(cluster fa)
outreg2 using "Table 2", excel adjr2 se alpha(0.01,0.05,0.10) ///
		symbol(***,**,*) bdec(3) sdec(3) rdec(3) ctitle("MultipleSons-p1") append

		
	
********************************************************************************
*********************************  THE END *************************************
********************************************************************************
